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1 Abstract 

This paper presents a simple, one-dimensional model of a randomly advected pas- 
sive scalar. The model exhibits anomalous inertial range scaling for the structure 
functions constructed from scalar differences. The model provides a simple compu- 
tational test for recent ideas regarding closure and scaling for randomly advected 
passive scalars. Results suggest that high order structure function scaling depends 
on the largest velocity eddy size, and hence scaling exponents may be geometry- 
dependent and non-universal. 
PACSs: 47.27.GS; 05.40.+j 
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2 Introduction 

Passive scalars are 'tracer particles' which are advected by the flow of another fluid. 
A typical example might be ink mixed by flowing water. The mixing of a scalar field 
T in a velocity field u is governed by the advection-diffusion equation, 

(d t + u(x, t) ■ V) T(x, t) = DV 2 T(x, t) + /(x, t) (1) 

in which D is a diffusion constant and / is an external source of scalar. The velocity 
field is determined independently of the scalar field (hence the adjective passive) and 
(in this discussion) is incompressible (V • u = 0). The incompressible flow simply 
transports the scalar, without changing any of the scalar's values. The diffusion 
term smooths the resulting scalar differences, while the source term generates new 
scalar to compensate for this smoothing. 

The mixing of a scalar field T is characterized experimentally by its structure 
functions, defined as S , 2n (r) = ((T(x + r, t) — T(x, t)) 2n ) for positive integers n. 
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It has been established experimentally that passive scalar structure functions are 
scaling functions of r over a wide range of length scales when the advecting flow is 
turbulent. This range of length scales, known as the inertial range, extends from 
a macroscale (usually determined by the flow geometry) down to a dissipative cut- 
off scale. The lower end of the inertial scaling range is set either by the viscous 
scale of the velocity field (the scale where viscosity damps out the fluid motion) or 
by the scalar's dissipation scale (the scale where molecular diffusivity damps out 
scalar fluctuations). The ratio of these length scales, known as the Prandtl number, 
is a material property of the system. This work will focus exclusively on systems 
in which the molecular diffusivity imposes the lower bound on the inertial range 
(Prandtl number less than one). In such systems the size of the inertial range 
is measured by the Peclet number, defined as the ratio of the macroscale to the 
dissipation scale. 

The inertial range scaling exponents p 2n , defined by S 2n (r) ~ r P2n , are mea- 
sured experimentally. The simplest possible prediction is based on the belief that 
the scalar difference between two points separated by a distance r will typically 
have a magnitude AT(r) which scales with r, AT(r) ~ r@ (for some j3). This im- 
plies S 2n {r) ~ r 2n /3. Consequently, the scaling exponents satisfy p 2n = np 2 , which 
is known as 'regular scaling.' Experimentally, the exponents are seen to exhibit 
'anomalous scaling,' p 2n < np 2 . Such scaling indicates that there is no typical 
magnitude for AT(r), and this wide variation in magnitudes is known as multi- 
scaling. Experimentally, measured scaling indices for passive scalars in turbulent 
flows are quite anomalous. |l|] Until recently, relatively little attention was paid to 
this observed anomalous scaling in the passive scalar field. This is perhaps because 
the phenomena had often been attributed to anomalous scaling of the underlying 
velocity field. 



However, a model proposed by Robert Kraichnan |13[ suggests that a hypo- 
thetical regular- scaling velocity field would still generate anomalous scaling in the 
scalar structure functions. In addition, the scaling exponents would be universal, 
independent of details such as the geometry of the flow. In this model, the pas- 
sive scalar is advected by a stochastic velocity field with regular scaling in space 
but with an extremely short correlation time. It is possible to form equations for 
the structure functions in this model, provided a closure assumption is made. The 



proposed closure [|13[] |L4j leads to structure functions which scale with r, but with 
scaling exponents p 2n which grow like yfn rather than n. Others, however, have ap- 
plied different techniques to the same model (avoiding a closure ansatz) and claim 
to have found different results. |@||A review of the disagreement is given 



by Shraiman and Siggia. [16 



To study this closure ansatz and the passive scalar structure functions numeri- 



cally, I use a one- dimensional model []10| in which the structure functions S 2n obey 
statistical equations with the same closure problem as in Kraichnan's model. The 
model is simple enough to allow direct numerical study of the structure function 
scaling (over 2 orders of magnitude in the separation r) and the proposed closure 
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ansatz. Mixing occurs through a random mapping function chosen to mimic an 
incompressible flow and induce scaling in the structure functions. Because it is re- 
stricted to one dimension, the model is best viewed as representing scalar mixing in 
a turbulent pipe flow. 

In this model the second order structure function can be determined analytically. 
For large Peclet number it scales with r with a scaling exponent p 2 determined by 
the advective mapping function. The higher order structure functions S , 2 n (r) also 
scale (numerically), with exponents p 2n which appear to approach a constant value 
at large n. 

The closure ansatz proposed by Kraichnan is seen to fail in numerical simulations 
of this model. The source of this failure can be traced to the finite size of the largest 
mixing events. This upper length scale is analogous to the largest eddy size in a 
fluid flow (in a pipe this size is determined by the pipe diameter). Very large scalar 
differences cannot be produced by a single mixing event, but instead arise from the 
combined action of several events. For this reason large scalar differences are expo- 
nentially improbable. Deviations from the closure ansatz are observed only for these 
large scalar differences. However, the behavior of the large n structure functions is 
determined by these rare events. In this model the scaling exponents approach a 
constant at large n, in contrast to the \fn behavior predicted for Kraichnan's model. 
In addition, the value of the constant might depend on the largest eddy size, so that 
the large n behavior of the exponents would be geometry dependent. 

Because these conclusions rest primarily on the fact that the largest mixing events 
are much smaller than the system size, one might expect them to be independent of 
the details of the small-scale mixing process. Mixing in a turbulent pipe flow should 
exhibit this separation of length scales, and there is some experimental evidence to 
support the idea that large scalar differences are exponentially unlikely in such a 
flow. || 



2.1 Kraichnan's Passive Scalar Model 

Robert Kraichnan first introduced the 'white-advected' passive scalar model in 1968 



T2I], and introduced the anomalous-scaling solution in a pair of recent papers. [T3 



14j] All of the principal results of this section can be found in the recent literature, 
although some of the notation has been altered here. 

In this model the passive scalar T(x, t) obeys the usual advection-diffusion equa- 
tion, 

{d t + u(x, t) ■ V) T(x, t) = DV 2 T(x, t) + /(x, t) (2) 

in which D is a diffusion constant. The external source / was not explicitly consid- 
ered by Kraichnan, but is necessary to maintain a state of statistical equilibrium. 
The incompressible velocity field (V • u = 0) is random with zero mean and is white 
in time. The second order velocity structure function is specified as 

([u(x + r, t) - u(x, t)) • -[u(x + r, t>) - u(x, t')]--) = CS(t - t>) (f-V (3) 



3 



where C is a constant and r] is a parameter of the model. This scaling in space 
holds until r approaches some large correlation length L v , representing the largest 
motions in the system. The scaling correlation in space is similar to real turbulent 
flows, while the very short range of temporal correlations is unphysical but desirable 
for technical reasons. The velocity statistics are chosen to be gaussian, so higher 
order velocity structure functions can be expressed in terms of this second order 
function. Consequently, the velocity structure functions exhibit regular scaling, 
and so any anomalous scaling found in the passive scalar structure functions must 
arise from the structure of the scalar advection-diflusion equation and not from the 
underlying velocity scaling. 

The source term is also a gaussian random variable with zero mean. Its correla- 
tions are specified by 

(f(x + r,t)f(x,t')) = X (r)6(t-t'). (4) 

The spatial function x(r) is assumed to be smooth, so that x( r ) — x(0) f° r r ^ L s , 
where L s represents the correlation length of the source field. 

Calculation of the structure functions proceeds by forming an equation for scalar 
differences at equal times, defined as A(x, y; t) = T(x, t) — T(y, t). One can take the 
ensemble average of this equation by integrating over the recent history of the scalar 
and using the defined correlation functions for u and /. The ^-correlation in time 
makes it possible to do the resulting integrals. Using the definition S 2n ('X- — y) = 
(A 2n (x, y)), the structure functions are found to obey 

(d t + C) S 2n = J 2n + F 2n . (5) 

In this equation, £ is the Richardson eddy-diffusivity operator, defined as 

CS 2n (r) = ^CL-^d^r^drS^r)) (6) 

where r = |x — y| and d is the dimension of space. The source term 

F 2n ee 2n(2n - l)S 2n _ 2 ( X (0) - x (r)) (7) 

was makes it possible to establish a state of statistical equilibrium (d t S 2n = 0). The 
dissipation function 

J 2n = 2n J D(A 2 "- 1 (V^ + Vj)A) (8) 

must be expressed in terms of S 2n before the system of equations for S 2n will be 
closed. 

The lowest non-trivial structure function, n = 1, can be determined exactly. The 
dissipation function J 2 can be evaluated by commuting derivatives with A and with 
the averaging process (...), so that 

J 2 = 2DV 2 r S 2 (r) - AD((VT) 2 }. (9) 
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The mean-squared scalar dissipation D((VT) 2 ) is constant, and can be determined 
by balancing against the scalar source (since the flux of scalar is conserved). At very 
large separations r ^> L s , the source correlation vanishes (x(r) — > 0) and S 2 must 
approach a constant. Consequently, the equation for S 2 implies 2_D((VT) 2 ) = x(0). 
The complete equation for S 2 can then be re-cast as 

Cr l - d d r {r d ' 1 (J~y d r S 2 {r)) + Dr l - d d r {r d - l d r S 2 {r)) = x(0) (10) 

assumming x(r) ~ x(0)- The two differential operators which act on 5*2 possess 
different scalings. The length scale at which the operators are equal is defined as 
the dissipation scale, r^: 

r d = L v {-) . (11) 



The upper limit of the inertial (scaling) range is set by the lesser of the two 
macroscales, L s and L v . In the inertial range (rj < r < min(L s , L v )), the sec- 
ond operator can be neglected, and the solution is approximately 

S2(r) K m$ fry (12) 



Cp 2 d \L V 

where p 2 = 2 — 77 is the second-order scaling exponent. Alternately, in the dissipation 
range (r -C r^), the first operator can be neglected, and the solution is approximately 

S 2 (r) ~ ^r 2 . (13) 
v ; 2Dd K J 

There are two possible approaches for determining the higher order scaling expo- 
nents. The first is to balance the eddy-diffusivity operator against the forcing term. 
If the 2n-th structure function has scaling exponent p 2n , then the corresponding 
eddy-diffusivity term scales with r with the scaling exponent p 2n — p 2 . The source 
term has scaling p 2n - 2 . If these terms are equal (i.e. if J 2n is negligible) then the 
scaling exponents are given by p 2n = np 2 , which is regular scaling. 



The other possibility, first proposed by Kraichnan fL3| , is to balance the eddy- 
diffusivity operator against the dissipation term J 2n . This basically assumes that 
the dominant scaling is set not by the source term, but by a homogeneous solution 
(or zero mode) of the structure function equations. Consequently, the scaling could 
be universal, independent of the details of the source term. 

Unfortunately, evaluation of the J 2n requires a closure ansatz stating how the 
laplacian of the scalar field at two distinct points, (V 2 , + V 2 )A(x, y) depends on 
the scalar difference A(x, y) between those points. The proposed closure ansatz is 
best understood in terms of the probability distribution function (PDF) for scalar 
differences P(A,r). This function gives the probability for seeing a particular value 
A for the scalar difference between two points separated by a distance r. The 
structure functions S 2n (r) are the moments of this PDF: 

S 2n (r) = J P(A,r)A 2n dA. (14) 
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The unknown closure information is expressed in the conditional probability 



H(A, x - y) = ((Vl + VJ) A(x, y) | A(x, y)> (15) 

which gives the ensemble-averaged value of (V^ + V^)A(x, y) subject to the con- 
straint that the scalar difference between the two points take on the particular value 
A. The dissipation functions can be expressed using H and the PDF for scalar 
differences P(A,r): 

J 2n {r) = 2nD J P(A, r)A 2n - x H{A, r)dA. (16) 



The closure ansatz proposed in |TJj is that H is approximately a linear function of 
A: 

H(A,r) ^ a(r)A (17) 

with some unknown slope a. This assumption can be viewed as a truncation of the 
Taylor series expansion of H(A) (by symmetry, the series has only odd terms). By 
substitution, the dissipation functions are J 2n = 2nDaS 2n . The function a can be 
determined by considering the n = 1 case. The result is that 

■hn = nJ 2 ^ (18) 

(recall J 2 is approximately constant in the inertial range). Hence both J 2n and CS 2n 
have scaling exponent p 2n — P2- Consequently, the scaling exponents are set by the 
coefficients of the terms. The result is that 

P2n = \\j^ndp 2 + (d- p 2 ) 2 - ~{d - p 2 ). (19) 



The asymptotic behavior is that p 2n — y/ndp 2 , which depends only on the spatial 
dimension d and the velocity scaling rj (through p 2 = 2 — 77). 



2.2 Competing Calculations 

The anomalous scaling for p 2n in the Kraichnan model depends on the closure ansatz, 
H(A) oc A. Unfortunately, it is not possible to consider higher order corrections to 
the Taylor series for H(A) because the resulting equations for S 2n cannot be solved. 
However, others have approached the model using techniques which avoid the need 
for a closure ansatz and, in certain limits, have found values of p 2n which disagree 
with Kraichnan's result. 

The competing approach to determining scaling exponents in this model focuses 
on the n-point correlation functions of the scalar (known as nth order moments), 
defined as 

M n ( Xl ,x 2 , ...x n ) = (T( Xl )T(x 2 )...T(x n )). (20) 
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The nth order structure function can be expressed in terms of the corresponding 
moment by allowing points to fuse together. The moments convey more information 
than the structure functions and can be viewed as the more fundamental objects, 
but moments are generally not measured experimentally 

The (^-correlation in time makes it possible to construct closed equations for the 
moments (known as Hopf equations). There is no closure problem associated with 
the dissipative terms. Unfortunately, the equation of the nth order moment is a 
partial differential equation in n\ vector separations, and it is beyond the reach of 
known mathematics to determine the solutions of such complicated equations. 

To determine scaling exponents, several groups have assumed an overall scaling 
for the nth order moment when all separations lie within the inertial range, and then 
attempted to calculate the scaling exponent in certain limits using perturbation 
theory. In each case it is assumed that the dominant contribution comes from a 
zero mode of the Hopf equation. These exponents differ from the structure function 
scaling exponents proposed by Kraichnan. The most readable description of this 
approach can be found in the paper by Chertkov, et. al. 

One of the first of these calculations was done in the limit p 2 — > 2. |J The result 
was that the 4th order scaling exponent was given by 

H - 2p 2 - ^ (21) 

(rj = 2 — p 2 <C 1). This indicates that the scaling is nearly regular in this limit, with 
a small correction for p 2 < 2. The solution proposed by Kraichnan (equation |T9| ) is 
not regular in this limit; instead it yields 

( 1 /T5 ; 1 A 1 / 2 + 3d \ 

p 4 = I + -d )--[ , =-l)n 22 

which is quite anomalous as rj —>■ (in three dimensions p^ — > 3, not 4). This 
disagreement is very substantial. Unfortunately, the huge difference in approach 
makes it difficult to ascertain the source of the disagreement. 

Another limit in which a descrepancy was found is the limit of large space di- 
mension [d — > oo). P| || In this case the scaling is again nearly regular, with a 
small correction in 1/d. The corresponding limit of equation is also regular, but 
the correction term differs from that obtained from the Hopf equation calculation. 

One proposed reconcilliation between these calculations is a generalization of 
Kraichnan's closure based on the 'fusion rules.' |5j The dissipation functions are 
given by a generalized form of equation 



■hn = nC n J 2 -^- (23) 

<->2 

with an unknown constant of proportionality C n . Values C n ^ 1 would change the 
result for the scaling exponents to 



P2n = -\J4nC n dp 2 + (d- p 2 ) 2 --(d-p 2 ). (24) 
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No precise form for C n has been proposed, but Ching |5| argues that C n oc n for 
large n, indicating p2 n oc n for large n. 

Numerical simulations of the Kraichnan model in two dimensions have been 



attempted by two groups. [14j [0] Although the size of the inertial range is limited, 
both groups claim to support Kraichnan's solution for intermediate values of pi- 
The limit p<i — > 2 has not been approached in either case. 



3 A Simple One-Dimensional Model 
3.1 Definition of the Model 

To address the debate, I propose a simple passive scalar model which reproduces 
the closure problem described above. The closure ansatz can be studied numerically 
on two levels. First, one can calculate the dissipation functions J2 n and measure 
the constants of proportionality to determine if C n ^ 1. Second, one can study the 
conditional probability H(A,r) directly to determine if it is a linear function of A. 

For numerical simulations, it is necessary to develop a model which is simple 
enough computationally to permit a large inertial range. Consequently, a one- 
dimensional scalar field is preferable. However, a one-dimensional incompressible 
velocity field would be quite dull, so we are forced to choose some other form of 
mixing that retains certain traits of an incompressible flow. The attribute preserved 
in this model is that the mixing simply re-arranges the scalar field, without altering 
any scalar values. Any advection of this type will necessarily be non-local. The par- 
ticular method of re-arrangement is based on a picture of a turbulent flow consisting 
of many swirling eddies, and is constructed to induce scaling in the passive scalar 
structure functions. Because the model is one-dimensional, it is best thought of as 
being analogous physically to a turbulent pipe flow. This model is loosely based on 
the 'linear eddy model' of Kerstein. JTT| 



The model is designed to produce structure functions which obey statistical 
equations with the same form as equation The dissipation functions J2n are 
identical, so that the closure problem is reproduced, but the eddy-diffusity operator 
C is altered (although it remains a scaling operator). 

The motivation for this particular model is developed as follows: Imagine that 
the one-dimensional scalar field is embedded in a plane. The advection consists of 
an eddy in that plane, centered on the scalar field, which rotates one-half turn. This 
maps the scalar field T{x) onto itself, according to the rule 

T(x) -> T(L - x) (25) 

where L is the eddy size (centered on x = \L). The eddy acts only on the region 
< x < L; elsewhere nothing happens. Applying one of these eddies in each time 
step r, with randomly chosen size and position, along with diffusion, gives a rule for 
advancing the state of the passive scalar by one time step, 

T(x, t) = T(x, 0) + V[T(x, 0)] + Drd 2 x T{x, 0) (26) 
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in which the advection operator V[T] is 

VfTfa^l — < r ^^ x ° ^ L — x) — T(x) if x a < x < x G + L, (27) 
[ otherwise. 

The model consists of applying this rule many times, with the size L and position x a 
chosen randomly at each step from appropriate probability distribution functions. 
All possible x a in the system have equal probability, but the eddy sizes L are chosen 
according to a scaling probability with scaling index y: 

P(L)dL = CL- y dL. (28) 

This generates scaling behavior within an inertial range determined by the smallest 
and largest possible values of L: L Q < L < L v . In the analogy to turbulent pipe 
flow, L v would play the role of the pipe diameter and L Q the role of a viscous cut- 
off scale (the size of the smallest motions). The passive scalar structure functions 
would be expected to exhibit inertial range scaling between L Q and L v , so long as 
the dissipation scale is not larger the L a . In practice, the constant Dt is chosen so 
that the dissipation scale is roughly equal to L (corresponding to a Prandtl number 
near unity), so that the ratio of length scales kf- plays the role of the Peclet number 
in the model. 

To maintain a state of statistical equilibrium, some forcing is required. In this 



model, the forcing is done by imposing an overall gradient g on the scalar, as in |T5 
A new variable 9 is then defined by the deviation from the gradient, 

9(x) = T(x)-gx (29) 

and structure functions are defined in terms of 9: 

S n (r) = ((9(x + r)-9(x)) n ). (30) 



3.2 Structure Function Equations 

Equations for the structure functions S n (r) can be constructed in the same way as 
in Section 1.2. The scalar difference A(x,y) = 9{x) — 9(y) obeys the equation 



A(x, y; t) = Dr{dl + d 2 ) A{x, y; 0) + *[A(z, y- 0)] 



(31) 



when time is advanced by one unit r. The action of the convective term $ on A 
depends on whether x or y (or both) lie within the eddy: 



V[A(x,y)) 



A(x,y) 

A(2x + L-x,y)+g(L + 2x c 



2x) 



x,y £ [x ,x + L\ 
x G [x Q , x Q + L\, 
y i [x , x + L] 

y e [x ,x + L] 

A(2x Q + L — x, 2x a + L — y) — 2g{x — y) x, y £ [x a , x + L). 



A(x, 2x Q + L - y) - g{L + 2x Q - 2y) 



(32) 
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Raising this equation to the nth power and taking the ensemble average gives 

S n {x-y-r) = (^[A(x,y;0)])+nDr(^- 1 [A(x ! y;0)](^ + ^)A(a;,y;0)) + ... (33) 

keeping only the lowest order term in Dr. In the diffusive term, \P[A] can be replaced 
with A, because only a small portion of the scalar field lies within the eddy at that 
time step (the eddy is another higher order correction). Then the equation becomes 

S n (x -y,r) = (* n [A{x, y- 0)]> + J n (x, y- 0) (34) 

where 

J n (x, y) = nDr(A(x, y) n -\d 2 x + d 2 y ) A(x, y)) (35) 

as in equation |8| and jr3|. The quantity (^ n [A(x,y; 0)]) can be computed from the 
definition of \I> by taking the ensemble average over the eddy probabilities for x Q and 
L. Defining the difference variable r = x — y, the result divides into three regions: a 
diffusive interval (r < L ), the inertial range (L Q < r < L v ), and a large scale region 
(r > L v ). The structure functions obey the equation 

A(S n (r, r) - S n (r, 0)) + C[S n (r, 0)] = AJ n (r, 0) + F n (r, 0) (36) 

where A is the system size. In statistical equilibrium S n (r,r) = 5 n (r, 0), and this 
equation is identical in form to equation [5| The eddy-diffusivity operator £ and the 
source F n (r) differ according to the value of r. For r < L a , they are 

C[S(r)] ee - [ Lv P(L)dL{ [ L+r dzS(z) - (r + L)S(r) - (r - L)S(-r)} (37) 

J h J L—r 

n ( m \ f L v r L+r 

F n (r) ee u ' \ i / P{L)dL{ / dzS n ^ m (z)(z - r) m 

~y m\(n — m)\ Jl Jl-v 

+ (L-r)5 n _ m (-r)(-2r) m }. (38) 
In the inertial range, L Q < r < L v , they are 

£[S{r)\ ee - T P{L)dL{ [ r+L dzS(z) - 2LS{r)} 
- [ V P(L)dL{ [ L+r dzS(z) - (r + L)S(r) - (r - L)S(-r)} (39) 

Jr J L—r 

F n(r) = E „ ' m { / / dzS n „ m (z)(z - r) m 

^ l m\[n — m)\ Jl Jt-l 

+ f Lv P(L)dL{ [ L+T dzS n _ m {z){z - r) m + {L- r)S' n _ m (-r)(-2r) m }}. (40) 

Jr J L—r 

In the large scale region, r > L v , they are 

C[S(r)] = - P{L)dL{ f +L dzS(z) - 2LS(r)} (41) 

J L Jr—L 

F n(r) = E u M / P{L)dL{ dzS n - m (z)(z-rr}. (42) 
r ~ l m\[n — my. J l Jr-L 
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3.3 Solution for S 2 (r) 

The first order structure function Si vanishes, so the lowest non-trivial structure 
function is S 2 . It can be determined approximately in both the inertial and dissipa- 
tive ranges. As a first step, it is necessary to re-write J 2 by commuting derivatives 
and using spatial homogeneity, as in the Kraichnan model: 

J 2 = 2Dr(A{d 2 x + d 2 y )A) = 2Drd 2 r S 2 - 4Dr((«9^) 2 > (43) 

where Dr({d x 9) 2 ) is the mean-square dissipation of the scalar (a constant). This 
constant can be evaluated by looking at the large scale region (where S 2 must ap- 
proach a constant) and balancing it with the forcing term F 2 = |g 2 (L 3 ), so that 

me?) - (44) 

This result serves as a check on the accuracy of numerical simulations. 

The solution far into the dissipative region (r <C L Q ) can then be evaluated by 
neglecting both C and F (setting J 2 = 0): 

S 2 (r) = ({d x ef)r 2 = lj^(9r) 2 . (45) 

In the inertial range, the approximate solution is found by balancing the convec- 
tive term C against the dissipation term J 2 = —4DT((d x 9) 2 ): 

C[S 2 (r)} = - 2 -(L 3 )g 2 . (46) 
By assuming a scaling solution, S 2 (r) = A 2 r p2 , the convective term becomes 

C[A 2 r^} = -^I(p 2 , y)r 2+ ^-y (47) 

in the limits ^ — > and ► 0. The definite integral I(p, y) is defined as 

I(p, y) = —{ / -((1 + zY+? - (1 - zf+P - 2(1 + p)z) 

1 + p Jo zy 

+ 1 ^^((1 + z) 1+P " (1 " *) 1+P - 2(1 + PV +P )} (48) 

and can be evaluated numerically. The first integral in I diverges at the lower limit 
for y > 3, setting an upper limit on the scaling range. (For y > 3 the lower limit 
must be replaced by — , and the ultraviolet region determines the solution.) Since 
£[5(r)] balances a constant, the scaling is fixed at p 2 = y — 2 (which only makes 
sense for y > 2) and the solution for L r L v is 

S 2{r) = 1. (£^f ■ ff r (49) 

where (L 3 ) has been evaluated explicitly. 

These solutions for ^2 compare very well with the numerical results in both the 
inertial and dissipative ranges (see Figure 0). 
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3.4 Higher Order Structure Functions 



Even-Order Structure Functions 

In the inertial range, the scaling of higher even-order structure functions is deter- 
mined by balancing C[S 2n ] against J 2n , as in |13| : 

C[S 2n (r)\ J 2n (r). (50) 

The closure assumption for J 2n is that 

S 

J2n{r) = nC n .l 2 -^-. (51) 

^2 

This generalized form with C n ^ 1 was suggested in 0, while C n = 1 gives the 
original Kraichnan ansatz of ||13|| . Assuming scaling functions for S 2n , the scaling 



exponents p 2n are determined by the coefficients in equation 

I(p2n, y) = nC n I(p 2 , y) (52) 



This equation for the exponents is analogous to equation |24| in the Kraichnan model. 
This result can be used (in principle) to evaluate the coefficients of proportionality 
C n given numerically measured p 2n . Unfortunately, the sensitivity of I to p limits 
the accuracy of this method. However, the numerical results for the exponents and 
C n (computed from equation |51]) are consistent with this result, indicating that the 
scaling is indeed set by balancing the convective term C against the dissipation term 

J2n- 



Odd-Order Structure Functions 

Odd order structure functions exhibit scaling only in the dissipative range. Because 
the eddy-diffusivity operator C differs depending on whether it operates on an odd 
or even function of r, the inertial range scaling solution exists only for even order 
structure functions. Odd order functions are positive scaling functions in the dissi- 
pation range (small r), but pass through zero in the inertial range and then approach 
zero from below at large r. This behavior of the odd structure functions is peculiar 
to this particular model and differs from physical passive scalars. 

3.5 Dissipation Range Scaling 

At small enough length scales the dissipative terms determine the scaling. In this 
(dissipative) range the J n term dominates the solution. It is convenient to re-write 
it as: 

J n {r) = 2Drd 2 r S n {r) - n{n - l)Dr(A n - 2 [(d x 9) 2 + (d y 9) 2 ]). (53) 

Balancing the two parts of J n against each other will lead to a solution if the con- 
ditional probability 

G(A,x-y) = (l(d x 6) 2 + (d y 6) 2 ]\A) (54) 
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can be calculated. In a purely dissipative system, G can be determined from a Taylor 
expansion of A at small separations: A = r(d x 9), which implies G = However, 
the presence of convection competes with the smoothing effects of diffusion and 
generates higher terms in the Taylor expansion of A. Given this fact, a reasonable 
closure approximation (which is supported numerically) is 

bA cA 2 

G(A,r) =a + — + — (55) 

where a, b, and c are constants. The assumed r dependence is necessary to generate 
a solution with regular scaling. 

For a regular scaling solution, S n (r) = A n r n , the unknown constants can be 
expressed in terms of A 2 (calculated above) and A 3 (unknown): 

A n = A n _ 2 A 2 + 4^A»-i. (56) 

A 2 

Hence all higher order structure functions are expressible in terms of S 2 and S3. 
Unfortunately, no analytic solution for S3 has been found. However, in the absence 
of any spatial assymmetry (which is generated by the gradient forcing used in this 
model) the odd order structure functions would vanish and the solution would be 

S 2n (r) = S£(r) (57) 
which has regular scaling but non-gaussian statistics. 



4 Numerical Results 
4.1 Simulation of the Model 

A computational model was created by discretizing the dynamical equation for T(x). 
The advective term is handled straightforwardly, since it is just a re-arrangement of 
the values of T(x). Diffusion was done subsequent to the advection process in each 
time step, using a second-order implicit finite differences scheme. An overall gradient 
was applied to generate forcing, and periodic boundary conditions for the fluctuating 
field 6{x) were used. The system was evolved to a state of statistical equilibrium 
before any averaging computations were done. Equilibrium was indicated by the 
establishment of stable (analytically known) values of ((d x 9) 2 ) and S 2 (r). Structure 
functions were computed by taking space-time averages over the entire simulation, 
implicitly assuming an ergodic system. 

Simulations were conducted in a system with 50000 grid points and a grid spacing 
of dx = 0.1. The smallest eddy size was L Q = 2, permitting a significant dissipative 
interval. The largest eddy size was either L v = 200 or L v = 500, allowing for inertial 
range scaling over two orders of magnitude. The imposed gradient was g = 0.01, 
and the diffusion constant Dt was of order 10" 5 (varying slightly with y to set the 
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Scalar Profile 
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A 'snapshot' of the scalar profile, revealing the overall gradient g = 0.01 along with 
significant fluctuations. The profile represents 6000 grid points (dx = 0.1). 

Figure 1: Typical Scalar Profile 



dissipation scale approximately equal to L Q ). The eddy scaling exponent was varied 
between y = 2.1 and y = 2.8 for L v = 200, and between y = 2.1 and y = 2.4 
for L v = 500. The number of time steps varied between 3 • 10 7 and 10 8 . A large 
number of time steps is needed to suitably average over the eddy size probability 
distribution P{L), especially for larger L v and y. In addition, several simulations 
were conducted with large diffusion constants (Dr = 0.02) to study the dissipative 
range solution. A typical 'snapshot' of the scalar profile T(x) is shown if Figure [I]. 
The profile reveals structure on a wide range of length scales, including significant 
'flat' regions, and appears similar to experimentally observed profiles. 

4.2 Numerically Computed Structure Functions 

Second order structure functions S2 have been computed analytically and can be 
compared directly with simulation results. Figure shows 5*2 as a function of r 
on a log-log plot (base e) for L v = 200 and four values of y. The smooth curves 
show the simulation data (50 data points each). The analytic results are shown as 
two straight lines (equation ^5] for the dissipative region, and equation [|9] for the 
inertial range). There are no free parameters. The agreement is good (within 5% in 
the heart of the inertial range) except at the inertial range boundaries, where the 
analytic calculation has no validity. Interestingly, the boundary terms apparently 
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Second Order Structure Functions 




Second order structure functions S2 (r) on a log-log plot, for L v = 200 and four values 
of the scaling parameter y. Smooth curves are numerical simulation data; straight lines 
represent the analytic scaling solution in two regimes (inertial and dissipative ranges). 



Figure 2: Second Order Structure Functions 



grow in importance as y increases, and this encroachment reduces the effective size 
of the inertial range. 

Higher even-order structure function also exhibit inertial range scaling. Figure 
|3] shows the even order scaling indices p 2n as a function of n for five values of 
y. The error bars represent the range of observed values over several simulations 
with different initial conditions. The scaling indices are independent of the upper 
length scale L v (within the error). The second order indices p2 lie within 3% of 
the theoretical values y — 2. The deviation from regular scaling (p2 n = npi) is quite 
pronounced. For larger y, the scaling exponents appear to approach a constant value 
(dependent on y) as n increases. 
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Inertial range structure function scaling exponents p2 n as a function of n for five values of 
y. These simulation data agree with the analytically known result p2 = y — 2 for n = 1. 
The data suggest that the exponents may approach a constant value as n — ► oo. 

Figure 3: Inertial Range Scaling Exponents 
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4.3 Probability Distribution Functions 

Figure ^ shows the probability distribution function (PDF) for A on a log-linear 
scale for several values of r from the y = 2.1, L v = 200 simulation. The core of the 
PDF, defined roughly by |A| < gL v (gL v = 2), is rounded. The slight asymmetry 
between positive and negative A is due to the imposed gradient. Figure |5] shows the 
PDFs for the same values of r from the y = 2.4, L v = 200 simulation. The core of 
this PDF is much more sharply peaked. 
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Probability distribution function .P(A) of scalar differences A for three separations r, from 
the y = 2.1, gL v = 2 simulation. The tails (|A| > gL v ) of the distribution appear to be 
exponential, with a slope which is independent of r. 

Figure 4: Scalar Difference PDF, y — 2.1 



All the numerically generated PDFs exhibit exponential tails for |A| ^> gL v : 

P(A,r) A(r)e" c|A| (58) 

where c is independent of the separation r. Exponential tails have also been de- 
rived in another model in a particular limit. J| This form of the PDF suggests 
that, for large n, the structure function scaling exponents p2 n approach a constant 
independent of n. In this case the structure functions obey 



lim 



2/i 



2n(2n - 1) 



5 



(59) 



2n-2 
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PDF of Scalar Differences 




Scalar Difference (Delta) 

Probability distribution function -P(A) of scalar differences A for three separations r, from 
the y = 2.4, gL v = 2 simulation. The tails (|A| > gL v ) of the distribution appear to be 
exponential, with a slope which is independent of r. 

Figure 5: Scalar Difference PDF, y = 2.4 



in the limit of large n. The PDF's and the structure functions can both be used to 
independently estimate c, and the results are shown in Figure 3.9 (as a function of 
p 2 ). The two sets of data represent the two different values of L v , L v = 200 and 
L v = 500. 

The exponential tails can be understood in terms of a random walk of fluid 
elements. To generate a particular temperature difference A, fluid elements initially 
separated by a distance of order A must be brought close together. Since the largest 
correlated motion in the system is of size L v , scalar differences larger in magnitude 
than gL v can only be generated by the uncorrelated action of several eddies. This 
behavior is essentially a random walk, and the multiplication of probabilities leads 
to the exponential PDF. 

To be more quantitative, consider the probability for the motion of a fluid element 
along a Lagrangian trajectory from position zero to position x. For simplicity assume 
that all eddies are the same size L. To move a distance x with (m — 1)L < x < mL 
(for integer m) will typically require the point to be moved by m eddies, each carrying 
it a distance of order L. However, the probability that the point will lie within a 
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particular eddy is 4, and these probabilities multiply: 



P(x ~ raV) 



A 



(60) 



The motion must occur quickly, in (order of magnitude) m steps, or else diffusion 
will cause the fluid element to equilibrate with its new environment. By assuming 
the two points involved in constructing A move independently, this probability can 
be converted to a PDF for A by using x ~ ^ regardless of the separation r. The 
result is 



~ e 



P(A) 



- 1 / f A 
C =jL ln {L 



-cA 



(61) 
(62) 



Assuming that each eddy is of size L v suggests estimates of c = 1.5 for L v = 200 and 
c = 0.5 for L v = 500. These values are shown as straight lines in Figure ^ and are 
in qualitative agreement with the data. Of course, not all eddies are of size L v ; in 
fact the number of large eddies decreases as pi increases. This results in the trend 
of increasing c shown in Figure |6|. 
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Slopes c of the tails of the PDF -P(A) measured from simulation data for L m = 200 and 
L m = 500. The straight lines represent the Lagrangian random- walk estimate of c. The 
dependence on p2 arises because large eddies become less frequent as p2 increases. 



Figure 6: Slopes of PDF Tails 
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Once two particles have come close together, they become subject to correlated 
motions and diffusion. This is presumably the source of the r dependence in the pre- 
factor A(r) of the exponential tail. Correlated motions have a weak influence on the 
PDF tail, and cannot alter its exponential character. Also, there is a slight asymme- 
try between A and —A in the prefactor; this is attributed to the fact that producing 
a negative A requires the two points to pass by each other (in one dimension), and 
during this time their motion is correlated. 



4.4 Kraichnan's Closure Ansatz 



To test the closure ansatz in |I3| and evaluate the constants of proportionality 
C n , the dissipation functions J 2n were computed directly from the simulation data. 
Rather than use the Laplacian, J 2n was re-written by commuting derivatives as 

J 2n (r) = 2Drd 2 S 2n (r) - 2n(2n - l)Dr(A 2n - 2 [(d x 6) 2 + (d y 6) 2 ]). (63) 

In the limit Dt — > 0, only the second term on the right side remains. In actual 
simulations, the first term makes a finite contribution which makes it more difficult 
to determine the inertial range scaling. So in practice only the second term (ap- 
proximated using finite differences) was used as a surrogate for J 2n in the inertial 
range. 

In the inertial range, the J 2n are scaling functions of r with scaling indices q 2n 
which, by the closure ansatz, ought to satisfy 

<?2n = P2n ~ p2- (64) 

Figure [7] shows that this part of the closure ansatz holds, by plotting q 2n vs. p 2n — P2 
for the L v = 500 simulations. 

The constants C n can be computed using C n = anc ^ avera gi n g over the 

inertial range (C n is approximately constant over the region averaged). The resulting 
values to not appear to depend on any parameters of the system other than y (a 
weak dependence on L v is suspected but not detectable). Figure |8| shows computed 
values for C n averaged over several simulations. The deviation from C n = 1 is small 



but significant, contradicting the ansatz of [|13j. The growing values of C n at small 
P2 suggests near-regular scaling in this limit, at least for moderate values of n. This 
is consistent with the scaling exponents in Figure |3|. The intermediate case p 2 = 0.4 
lies very close to the Kraichnan prediction of C n — 1. For larger p 2 , C n decreases 
with n, consistent with the scaling exponents approaching a constant value (as in 
Figure ||) . 

In ||14|| , the closure ansatz in [13| is derived from an assumption about a condi- 
tional probability, namely that 

ff(A, x — y) = ((fig + 6J) A(x, y)\A(x, y)) = (^r) A. (65) 
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Scaling Exponents Comparison 
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Inertial range scaling exponents q 2n for the dissipation functions J 2n compared to the 
structure function scaling exponents p 2n . Kraichnan's closure ansatz requires q 2n = P2n — 
P2, and the simulation data is consistent with this result. 

Figure 7: Scaling Exponent Comparison 



This conditional probability H(A, r) has be computed numerically in the simula- 
tions. Figure |9] shows H(A) as a function of A, normalized by 2 dtS 2 ' ^ or P 2 = ^ 
and both L v = 200 and L v = 500. Two values of r are shown for each simulation. 
The resulting averages lie very close to H(A) = A for small A, as assumed in Jl4| . 
However, they deviate from the straight line at |A| = 2gL v . This is approximately 
the point at which correlated motion becomes unimportant and the dynamics are 
controlled by the 'random walk' described on the previous section. The failure of 
Kraichnan's ansatz appears to be due to the existence of a finite upper size of the 
inertial range, L v , which is much smaller than the system size. Other values of p2 
exhibit a similar behavior. 

Upon closer examination, the large n behavior of the constants of proportionality 
C n can be related to the failure of the closure ansatz. Using the definition of the 
constants, 



C n = 



(66) 



nJ 2 S 2n 

the C n can be related to the closure ansatz through the definitions of S 2n and J 2n 
These definitions are 

(67) 



S 2n = J P(A)A 2n dA 
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Closure Constants 
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Closure constants of proportionality C n , defined by J^n = nCnJz-g^, from simulation data 
for several y. Kraichnan's closure ansatz requires C n = 1, which is inconsistent with the 
data. For larger y, C n decreases with n, suggesting scaling exponents might approach 
a constant value at large n. For smaller y, C n increases with n, suggesting near-regular 
scaling for moderate values of n. Both trends are consistent with the scaling exponents 
shown in Figure |3|. 

Figure 8: Closure Constants of Proportionality 
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Conditional Probability Closure Ansatz 
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Conditional probability H(A,x — y) = ((d^ + dy)A(x,y)\A(x,y)) as a function of the 
temperature difference A for two values of L v and two separations r = x — y. Data from 
the y = 2.1 simulations. H is normalized so that H = A is the Kraichnan ansatz. The 
ansatz is in good agreement with the data for small values of A, but the data deviates 
for |A| > 2gL v {g = 0.01). This is approximately the largest value of A which can be 
generated by a single large eddy; larger values are exponentially unlikely (see Figure 

Figure 9: Conditional Probability Closure Ansatz 
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and 

J 2n = nDr J P(A)A 2n ~ 1 H(A)dA (68) 

where P{A) is the probability distribution function for scalar differences. 
If one defines the 'error' in Kraichnan's closure as 5H 

8H(A) = A- ^^H(A) (69) 

then the closure constants can be written as 

C = ^ (70) 

where the functions Q 2n are defined as 

^ /P(A)A 2 "-^A 

representing the contribution of the 'error' SH to C n . 

Numerically, 8H is aproximately zero for small | A | , but appears to grow linearly 
with A for |A| > 2gL v . Asymptotically, H(A) might grow as fast as linearly with 
A (but with a slope less than one). An approximate large |A| form for H consistent 
with the simulation results is 

H ^- a 23k A < 72) 

where < a < 1. Then 5H oc (1 — a)A for |A| > 2gL v . Assuming this form for 
5H, the function Q 2 should be less than one. However, as n — > oo, Q 2n — > 1 — a. If 
the numerically observed trend for H(A) persists as |A| — >• oo, this implies 

lim C n = (73) 

which is a constant. It would be consistent with the numerics to have a = and 
hence 

lim C n = 0. (74) 

n— >oo 

It is also possible that a is non-zero and that l _ a q 2 depends on p 2 . This might 
explain the possibly different asymptotic values for C n for different values of p 2 seen 
in Figure |[ 

Unfortunately, in the Kraichnan model the truly interesting quantity is nC n . For 
a = 0, there are two possible assymptotic behaviors. If nC n is bounded at very large 
n, the scaling exponents p 2n would approach a constant as n — ► oo (as suggested for 
the eddy model). If instead nC n is unbounded (but growing no faster than n), then 
p 2n would also be unbounded, but growing no faster than y/n. The possibility that 
nC n — ^ is excluded by the Holder inequalities, since it would imply p 2n decreasing 
with n. If instead a > 0, nC n must grow linearly with n and the scaling exponents 
ultimately grow like y/n, although the precise values may differ from the C n = 1 
prediction. 
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4.5 Dissipative Range Closure Ansatz 



The conditional probability G(A) needed for closure in the dissipative range is shown 
in Figure |10| for several values of r (from the simulation with Dt = 0.02). The 
conditional probabilities are well approximated by the parabola 



G(A,r) 



bA cA 2 
a H 1 7T- 



(75) 



with a = 4 ■ 10 7 , b = 2 • 10 3 , and c = 1.99. These values are consistent with the 
independently known values of A 2 and A3. 

Dissipation Closure Ansatz 
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Conditional probability G(A,x — y) = ([(d x 9) 2 + (d y 9) 2 ]\A(x,y)) as a function of the 
temperature difference A for several values of r = x — y in the dissipative subrange. 
Data from the diffusion-dominated simulation. The parabolic form is consistant with the 
small-scale closure ansatz proposed in the previous chapter. 



Figure 10: Dissipation Range Closure Ansatz 



The structure functions are given by S n (r) = A n r n , and the constants A n are 
shown in Figure [11]. In addition, the calculated values A n = A n _ 2 A2 + j^A n _i from 
equation KBI are also shown, using the analytically known value of A 2 and the value 
of A3 determined from fitting 5*3 (r). The deviations from the analytic solution at 
large n are possibly due to the neglected terms F n and £ in the equation for the 
structure functions. 
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Dissipative Range Solution 



J -50 




Structure function coefficients Ai n as a function of n for the dissipative range regular- 
scaling solution, S2n = A2 n r 2n . The data points were determined by fitting the simulation 
data. The straight line is the analytic solution, with the unknown parameter ^| deter- 
mined by fitting the numerical data for Ss(r). 



Figure 11: Dissipation Range Solution 
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5 Conclusions 



5.1 Extensions to the Kraichnan Model 

The original motivation for studying this 'eddy model' was to test the validity of 
the closure ansatz proposed in JL4|. In the model, the principal failure of the closure 
could be traced to the finite size of the largest eddy, L v . Specifically, the conditional 
probability H(A) is not a linear function of A, deviating at large values (|A| 3> gL v , 
where g is the imposed gradient). The deviation arises because the largest velocity 
eddy in the system (L v ) is much smaller than the system size, so that large scalar 
differences only occur as a result of the uncorrelated action of several eddies. 

This finite L v effect causes the constants of proportionality C n proposed in |5| 
to deviate from 1, growing with n for small p2 but shrinking with n at large p2 (in 
contradiction to the predictions of || and also the simulations of @). Other apparent 
consequences of this effect include exponential tails for scalar difference PDFs and 
scaling exponents p2 n which approach a constant value at large n. Because this 
effect depends only on the existence of a finite eddy size L v and not on the details 
of the smaller scale mixing, it is reasonable to expect that it might apply generically 
to other models of scalar mixing, including the one proposed by Kraichnan. 

The limit L v —> oo 

For small values of the scalar difference A (|A| < gL v ) the numerical simulations 
of the eddy model appear to support the conditional probability ansatz H(A) oc 
A. Low order structure functions should therefore have exponents p2 n which lie 
very close to the values predicted by Kraichnan. The fact that the errors become 
significant at large |A| means that the scaling exponents pm will exhibit deviations 
from the values predicted by this ansatz when one looks at sufficiently high order 
structure functions (large n). 

Mathematically, structure function scaling is defined for an infinite inertial range, 
oi L v — > oo. When considering the large n behavior of the scaling indices, one 
must take the limits L v — > oo and n — > oo. The order matters. If one considers 
small enough values of n (for any particular L v ) it will appear that Kraichnan is 
correct. Hence if one lets L v — > oo first, all values of n will satisfy his predictions. 
For moderate values of p2 the numerical simulations of the 'eddy model' support 
Krachnan's ansatz if the limits are taken in this order. However, for any finite L v 
one can find sufficiently large n for which the scaling exponents deviate from the 
Kraichnan values due to the finite L v effect. Hence if one lets n — > oo first, the 
Kraichnan solution will fail. It is my personal opinion that keeping L v finite while 
letting n — > oo is the more physical limit, since any real system will have a finite 
largest eddy size. 



27 



The limit p 2 -> 

There is some numerical evidence that regular scaling is approached in the limit 
p 2 — > in the eddy model. The constants of proportionality deviate significantly 
from C n = 1 even for the lowest values of n, and the lowest scaling exponents fall 
just below the regular scaling line p 2n = np 2 . In addition, the core of the scalar 
difference PDF becomes very rounded, perhaps gaussian, in this limit. The source 
of these effects is not understood, but it is very possible that they have nothing 
to do with the finite L v effect. If so, they may indicate a failure of the Kraichnan 
ansatz which would persist in the limit L v — > oo. However, it must be remembered 
that direct calculation of the conditional probability H(A) does not indicate any 
measurable discrepancy with the closure ansatz for small values of A in this limit. 

Comparison with other Calculations 

In most of the competing calculations for the exponents of the Kraichnan model 
Hi HI Hi) the gradient forcing is replaced by a random source term for the passive 
scalar. In this formulation there are two upper length scales in the problem. They 
are the largest correlation length of the velocity field, L v , and the largest correlation 
length of the source field, L s . Whether or not one would expect to see deviations 
from Kraichnan's ansatz of the type observed in the eddy model depends on which 
of these two length scales is larger. 

First, consider L s ^> L v . The source of the scalar is smooth over length scales 
much larger than the largest motions of the velocity field. In this case very large 
scalar differences can only be generated by uncorrelated transport over scales much 
larger than L v . This is analagous to the transport against the gradient by random 
events in the eddy model. Hence in this case one might expect to see scalar difference 
PDFs with exponential tails, and deviations from Kraichnan's ansatz like those seen 
in the eddy model. 

The opposite case is L s C L„. Then the largest 'globs' of scalar (of size L s ) can 
be mixed by coherent motions of the velocity field. The sort of transport described 
above would not happen, and the Kraichnan ansatz could be correct in this case. 
However, this is the order of limits taken in P H, and these authors claim to 
find deviations from the Kraichnan result. If they are correct, the deviations come 
from a different source than the finite L v effect, so the eddy model simulations don't 
necessarily support their claims. 

A recent paper by Chertkov, et. ai, 0, has more in common with the eddy 
model. It considers a scalar advected by a one-dimensional compressible velocity 
field with correlations analogous to the Kraichnan model. By assuming that the tails 
of the scalar difference PDF are dominated by the most rapidly stretched Lagrangian 
trajectories, they conclude that the tails must be exponential and that the scaling 
exponents p 2n must be independent of n. The calculations were done in the limit 
Pi —>■ 0. Although the level of mathematical sophistication is quite different, the 
basic assumption is identical to the one used to understand the PDF tails in the 



28 



eddy model, and the conclusions are also similar. 
5.2 Extensions to Pipe Flow 

As suggested earlier, one can draw physical analogies between the eddy model and 
scalar mixing in turbulent pipe flow. In particular, the separation of length scales 
(between the largest eddy size and the system size) required for exponential tails 
in the scalar difference PDF can be achieved in pipe flow. The finite diameter 
of the pipe limits velocity eddies to this size, while the pipe itself can be much 
longer (analogous to this one-dimensional model). The mixing in the eddy model 
at small scales is quite different from real flows, but the likelihood of transport 
across distances larger than L v (analogous to the pipe diameter) could be similar 
(statistically) to real pipes. This suggests that, for large n, the scaling exponents 
in pipe flows would approach a constant value. Further, since the slope of the PDF 
tails depends on the pipe diameter (gL v ), this value might be geometry-dependent. 

Recent experiments on pipe flows |J have revealed PDF's of scalar values with 
exponential tails. No attempt has been made to study scalar differences, but the 
arguments of the previous chapter suggest that exponential tails in the scalar value 
PDF would imply exponential tails in the scalar difference PDF. If true, this would 
force structure function scaling exponents in turbulent pipe flow to approach a 
constant (non-universal) value at large n. 

Flows in other geometries (such as a box, for example), would not be expected 
to display this effect. If the largest velocity eddy is comparable to the system size, 
large-scale mixing can occur under the correlated motion of a single large eddy, and 
the mechanism for generating exponential tails in the PDF would not exist. Hence 
there is no reason to expect that scaling exponents pm would approach a constant at 
large n in such geometries. The structure function scaling exponents of the passive 
scalar might therefore be geometry-dependent, with pipes exhibiting a different large 
n asymptotic behavior than more open flow geometries. 
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